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Antisymmetrized molecular dynamics with quantum branching is generalized so as to allow finite 
time duration of the unrestricted coherent mean field propagation which is followed by the decoher- 
ence into wave packets. In this new model, the wave packet shrinking by the mean field propagation 
is respected as well as the diffusion, so that it predicts a one-body dynamics similar to that in mean 
field models. The shrinking effect is expected to change the diffusion property of nucleons in nuclear 
matter and the global one-body dynamics. The central 129 Xe + Sn collisions at 50 MeV/nucleon 
are calculated by the models with and without shrinking, and it is shown that the inclusion of the 
wave packet shrinking has a large effect on the multifragmentation in a big expanding system with 
a moderate expansion velocity. 

PACS numbers: 24.10.Cn, 25.70.Pq, 02.50.Ey, 02.70.Ns 



I. INTRODUCTION 

In order to describe heavy ion reactions as the dynam- 
ics of many-nucleon systems, two different kinds of mi- 
croscopic approaches have been proposed and applied. 
One is the molecular dynamics models jl], [|, [| % and 
the other is the mean field (TDHF-like) models g, §. 
An advantage of the mean field models is that they do 
not put any restriction on the one-body motion, while 
their disadvantage is that they cannot properly describe 
the cluster formation because of the lack of many-body 
cluster correlations. On the other hand, usual molecu- 
lar dynamics models assume a fixed Gaussian shape for 
single particle wave functions. This is an efficient way 
to describe the cluster correlation even by using a simple 
product wave function with or without antisymmetriza- 
tion. However, in another sense, the use of localized wave 
packets can be a regression because the one-body dynam- 
ics is not as precisely described as in mean field models. 

An unified understanding is desired on the question 
whether the single particle wave functions should be un- 
restricted or localized. Unless we can solve the dynam- 
ics keeping the full order of the many-body correlations, 
it is essential for a reasonable model to introduce the 
fluctuations which bring the system into many quantum 
branches each of which corresponds to one of the reaction 
channels or the configurations of clusterization. In mean 
field models, it has been proposed to introduce fluctua- 
tions in the one-body distribution function |?|, ^] , though 
it should be a difficult problem to determine the corre- 
lations among an infinite number of the degrees of free- 
dom of fluctuations. In this viewpoint, the philosophy 
of molecular dynamics is to introduce a special kind of 
fluctuation by stochastically localizing the single parti- 
cle wave functions, which is essential for the cluster pro- 
duction. The mean field equation should be interpreted 
as giving the short-time evolution of the one-body dis- 



tribution averaged over the stochastic branches. Based 
on this idea, the antisymmetrized molecular dynamics 
(AMD) has been extended in Refs. |, by incorporat- 
ing the wave packet diffusion effect in the mean field as 
a source of the fluctuation to the wave packet centroid 
which causes the quantum branching on the many-body 
level. Thus the single particle wave functions are local- 
ized in each branch, which makes cluster formation pos- 
sible, while the single particle motion is not restricted 
for the averaged value over the branches. The stochas- 
tic treatment of the dynamics of the wave packet width 
is the essential point of our approach. There is another 
approach j|, [o, |l2[ [r| which treats the width parame- 
ters as time-dependent variables in molecular dynamics, 
though it has turned out that the deterministic dynam- 
ics of the width variables cannot explain the evolution 
of the density fluctuation and the multiple cluster for- 
mation Jll|, [[2). Ohnishi and Randrup JlJ, |l5) have in- 
troduced quantum fluctuation into wave packet molecu- 
lar dynamics based on their statistical consideration and 
have shown its importance in cluster formation. However 
the dynamical origin of the their fluctuation has not been 
made clear. 

An unsatisfactory point of the improvement in Refs. 
H was that the stochastic fluctuation to the wave 
packet centroids can diffuse the distribution but cannot 
shrink the distribution. Note that a wave packet in the 
mean field normally diffuses in three directions in phase 
space and shrinks in the other three directions. Because 
of this difficulty, for example, the improved AMD could 
not be directly applied to an isolated nucleon, and there- 
fore the diffusion was switched off for isolated nucleons, 
which introduces an ambiguity to the model. To solve 
this kind of problems, we absolutely need a consistent 
understanding of both the mean field description and the 
molecular dynamics description. 

The first purpose of the present work is to construct 
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a general framework which contains both molecular dy- 
namics models and mean field models as specific cases. 
In this framework, the time evolution of a many-body 
system is given by the coherent mean field propagation 
and the decoherence of single particle states into wave 
packets. It will have two physically essential ingredi- 
ents, t and r m f, which define when "decoherence" and 
"mean field branching" take place, respectively. The 
choice (r, r m f) = (oo, oo) corresponds to a mean field 
model, while the choice (r, r m f) = (0, 0) corresponds to 
the version of AMD of Refs. || [ftj (called AMD /D in this 
paper) . Based on this general framework, we introduce a 
new model AMD/DS as the case of (r, r m f) = (large, 0), 
with which we can respect not only the diffusion but also 
the shrinking of the phase space distribution predicted 
by the coherent mean field propagation. 

The second purpose of the present work is to demon- 
strate the effect of the wave packet shrinking. The dif- 
ference between AMD/D and AMD/DS is expected to 
result in the different diffusion properties of nucleons in 
nuclear matter and the different global one-body dynam- 
ics. We perform the AMD/D and AMD/DS calculations 
for the central 129 Xe + Sn collisions at 50 MeV/nucleon 
and compare the results. The velocity of the expansion 
strongly depends on the model. The different expansion 
velocity results in the different cluster size distribution. 
It is shown, by comparison with INDRA experimental 
data, that AMD/D had problems of the overestimation 
of Z = 4, 5, 6 clusters and the underestimation of Z > 15 
clusters and these problems are solved by AMD/DS. 

This paper is organized as follows. In Sec. II, the for- 
mulation is presented. The physical principle of the gen- 
eral framework is giv en in II A, and then the concrete 
formulae are given in [IB which includes all the details. 
In |II C| , we introduce specific models such as AMD/D 
and AMD/DS, and give simple examples to show how 
our formulation works for AMD/D and AMD/DS. In 
Sec. Ill, the results of the calculations with AMD/D 
and AMD/DS are compared to each other and to the 
INDRA data for central 129 Xe + Sn collisions at 50 
MeV/nucleon, so as to demonstrate the important effect 
of the wave packet shrinking in AMD/DS in multifrag- 
mentation. Section IV is devoted to a summary. 



II. 



MEAN FIELD PROPAGATION FOLLOWED 
BY DECOHERENCE 

A. Principle 



First we give a general framework which includes both 
mean field models and molecular dynamics models as 
limiting cases of model parameters. Originally we have 
reached this framework in the course of an extension of 
AMD by the incorporation of good points of mean field 
models. However, the framework is given here from a 
more general point of view as an approximation of quan- 
tum many-body dynamics. 



In what follows, the two-nucleon collision effect is not 
shown explicitly for the brevity of presentation, but it 
is always considered in all the practical calculations as a 
stochastic process [[| ^6). 

In this model, we will use (without specifying yet how 
to use it) the AMD wave function which is a Slater de- 
terminant of Gaussian wave packets [9 , 
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exp{-^-^) 2 }xa 4 (i) 



The complex variables Z = {Zi 



(1) 

= 1, . . . , ^4} represent 
the centroids of the wave packets. We take the width 
parameter v = 0.16 fm~ 2 and the spin isospin states 
Xai — P 1\ P I) n Tj or n I- Because of the antisym- 
metrization, the variables Z do not necessarily have a di- 
rect physical meaning as nucleon positions and momenta. 
The AMD wave function |$(Z)) contains many quantum 
features in it and has been applied to the study of nu- 
clear structure problems with some extensions such as 
the parity and angular momentum projections |17| ]. 

The time-dependent many body wave function iNf'(i)) 
describing a complicated nuclear collision is, in an inter- 
mediate or final state, a superposition of a huge number 
of channels each of which corresponds to a different clus- 
terization configuration. We do not try to directly treat 
such a complicated many-body state nor we do 

not approximate \^(t)) by a single AMD wave function 
\§{Z)). We rather approximate the many-body density 
matrix y\> (t)) (t)\ by an ensemble of many AMD wave 
functions. 



m))m)\ 



Mz))(*(z)\ 

{*(Z)\*(Z)) 



w(Z, t)dZ. (2) 



A pure many-body state is approximated by a mixed 
state. This approximation itself will not be a serious 
drawback practically because the nuclear collision dy- 
namics is so complicated that one cannot observe full 
many-body correlations to distinguish a pure state from 
a mixed state. The benefit of this approximation is that 
the right hand side of Eq. (Q) still contains nontrivial 
many-body correlations required in multiple cluster for- 
mation, even though an ensemble of AMD wave func- 
tions is sufficiently simple to be tractable numerically. 
We should, of course, define a reasonable time evolution 
of the weight w(Z,t) or, alternatively stochastic trajec- 
tories of the variables Z(t). 

What is the physical idea when we use AMD wave 
functions \$(Z)) in Eq. (^)? When multiple cluster for- 
mation takes place, |\&(t)) is composed of a huge number 
of branches. The decomposition into branches should 
be done so that, in each branch, the one-body distribu- 
tion of each nucleon is localized in one of the clusters. 
Namely, a nucleon should not belong to many clusters 
at the same time, to avoid non-integer mass numbers of 
clusters. [The width of wave packets v has been chosen 
in such a way] In turn, if nucleons are localized in each 
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branch, the clusters are naturally bound due to the mean 
field among localized nucleons. This idea exists behind 
the molecular dynamics models which restrict each single 
particle wave function in each branch to a wave packet. 

The time evolution in our model is determined by two 
factors, the mean field propagation and the decomposi- 
tion into branches. At a time to, let us take one of the 
branches \§(Z)){$(Z)\ from Eq. (§). This is justified 
because the time evolution of a branch is independent of 
the others due to linear quantum mechanics. We consider 
the mean field propagation from t to t + r, 

|$(Z))($(Z)|-,|*(r,Z))<*(r,Z)|, (3) 

where |^(r, Z)) is the solution of the mean field equation 

l h±\*(t,Z))=H™(t)\*(t,Z)) (4) 

with the initial condition |*(0, Z)) = \§{Z)). The mean 
field Hamiltonian has a form 

H™{t) =£^)=E(|J + U(r h p 2 ; p(t))) , (5) 

i=l i=l ^ ' 

where the potential U depends on a one-body density 
matrix p(t). For the moment, we may assume that p(t) 
is the one-body density matrix for the state \^(t, Z)). We 
wish to emphasize that the single particle wave functions 
are not restricted to Gaussian packets in the mean field 
propagation and therefore |\P(t, Z)) is a general Slater 
determinant. Next, at the time to + r, the propagated 
state \^(t,Z))(^!(t,Z)\ is decomposed into AMD wave 
functions as 



|tt(r,Z))(tt(r, Z)\ 
(HZ)\$(Z)) 



(*(*)!*(*)) 



w(z,t; Z)dz, (6) 



with a suitable weight w(z,t; Z). A reasonable princi- 
ple to determine w(z, r; Z) would be first to choose a set 
of important (one-body) operators {O a } and then to re- 
quire that both sides of Eq. (^) give the same expectation 
values for all {O a }, 



($(z)|6 a |$(z)) 

(*(*)!*(*)> ' 



)(z, r; Z)dz = 



(*(r,Z)|6 a |*(T,Z)) 



($(Z)\*(Z)) 



(7) 



However, this prescription will not always work well 
because it is incompatible to the necessary condition 
w(z, r; Z) > 0. A practical choice of {O a } will be given 
bellow. In this way, in principle, Eq. (|3j) and Eq. (|6|) de- 
termine the time evolution from to to to + t. The next 
time evolution after to + t is obtained successively by 
applying the same model to each term of the right hand 
side of Eq. §). 

What is the physical meaning of r? Is Eq. (^) just a 
numerical approximation of a mean field model or does 
it have any physical meaning? Once one believes that 



the mean field propagation were always perfect, then the 
choice of a finite r would be unphysical. But the mean 
field propagation is not perfect at least in that it un- 
physically keeps the idempotency p 2 = p of the one-body 
density matrix. In many-body systems, even though the 
idempotency is satisfied at a time, it is not the case at 
a later time due to many-body correlations. The re- 
duced one-body density matrix will be an ensemble of 
idempotent density matrices. If we ignore the antisym- 
metrization for the simplicity of the discussion here, the 
one-body density matrix of a nucleon will be given like 



p = w'\ip'){ip'\ +w"\<p"){<p" 



(8) 



in which different components \tp'), \ L p"), ■ ■ ■ do not in- 
terfere in a simple way. The many-body state for this 
situation will be like 



|*) = cV> ® |*'> 



c'V'> 



I*'' 



(9) 



(§). In 



where the states of the rest of the system |^') 
. . . are orthogonal to one another, to have Eq. 
our model, the coherent mean field propagation 
\ip) ® \^f)) is assumed to be valid in the time duration 
from to to to + t, and "decoherence" is assumed to phys- 
ically take place into Gaussian packets at to +r as in Eq. 
d) and Eq. (§). [ "Decoherence" is a general concept of 
quantum mechanics in open systems |lq]. It means not 
only that Eq. (||) is satisfied at a specific time but also 
that different branches do not interfere at any later time.] 
In this way, the parameter r has a deep physical meaning 
as the coherence time. We do not try to determine r a 
priori in the present paper, but we assume that decoher- 
ence is a physical process as should be the case at least in 
multiple cluster formation. It should be noted that the 
physical state changes by decoherence, and therefore Eq. 
(0) should not be required for all the operators {O a }. We 
should require Eq. (M) for those operators {O a } which we 
can believe to remain after physical decoherence. 

Closely related to the decoherence of a nucleon, we 
should decide how this nucleon interacts with the rest of 
the system, namely, how this nucleon contributes to the 
mean field potential U. Once decoherence takes place at 
to+r so that different branches in Eq. do not interfere 
any longer, then the mean field should be calculated in 
each branch, independently of the other branches, by us- 
ing the corresponding wave packet (one of \ \<p"), ■ ■ ■ 
). We will call this change of the mean field "mean field 
branching" . The decomposition of the many body state 
in Eq. (|6|) is consistent to the mean field branching at 
to + t. Nevertheless, we may think of the other possibil- 
ity that the time scale of mean field branching, denoted 
by T m f, is shorter than r. The choice r m f < r can be 
reasonable in such a physical situation that, even before 
decoherence, the mean field approximation is applicable 
not to I*) of Eq. (||) but to each of the "pre- branches" in 
the right hand side of Eq. (||) even when the pre-branches 
have not been decohered yet (i.e., 1^'), I'J'"^ . . . are not 
orthogonal). This situation is possible because the true 
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time evolution is linear while the mean field approxima- 
tion is nonlinear. Therefore we regard r m f as a physical 
ingredient of the model which is not necessarily equal to 
r. 

The physical origin of decohcrence we have in mind 
here is the full oder of many-body correlations beyond 
the two-body correlation, which is especially important in 
multiple cluster formation. The choice of wave packets as 
decohered states is done for this purpose. Of course, the 
two-body collision effect destroys the idempotency of the 
one-body density matrix, but this effect is already taken 
into account by a more explicit way of the stochastic 
two-nucleon collision process || [T(| in all the realistic 
calculations. 



We wish to emphasize that much of the fermionic fea- 
ture remains even though we adopt the physical coordi- 
nate approximation [Eq. (pjj)]. This is because the value 
of W carries the fermionic information. For example, in 
the 6j4-dimensional space of W, there exists Pauli for- 
bidden region where W cannot take value for any choice 
of the original coordinate Z ||. Nevertheless, it depends 
on the purpose whether this approximation gives a good 
result. For example, it gives only a poor approximation 
for the evaluation of the Hamiltonian and the mean field 
potential. For the calculation of such quantities, we use 
the exact antisymmetrized wave function |$(.Z)) or use 
a better approximation scheme llQj. 



B. Formulation 

Now we start to give concrete formulae for the calcu- 
lation of the time evolution governed by the mean field 
propagation followed by the decoherence into stochastic 
branches. The explanation will be given by five steps. 

Impatient r eader s may be able to skip this part first 
moving to Sec. |EI Q, and later come back here if necessary. 



1. Physical coordinate approximation 

We adopt the "physical coordinate approximation" for 
an AMD wave function, 



MZ))(*(Z)\ 



\<p(w k ) Xak )(vCw k )x ak 



(10) 



where the spatial wave function of each nucleon k is given 
by a Gaussian packet 

(r|^(W fc ))^exp{-^(r-^) 2 }. (11) 

The centroids are the physical coordinates W = {Wfc} 
defined in Ref. § by 



V^ + ^P fe = £(v/Q) fcj Z;, (12) 



Qkj = BkjB jt ^, Bkj = e Zk ' Z:, S akaj . (13) 



In the phase space representation, the Wigner function 
for the nucleon k is given by a Gaussian packet 

6 



g(x;X k ) = 8exp 



(14) 



where we have introduced the 6-dimcnsional phase space 
coordinates 



(15) 
(16) 



x = {x a ; a = l,... ,6} = l^/vr 
X k = {X ka ; a = 1, . .. ,6} = \ V^Rfc 



2%sjv J 



2. Mean field propagation 

The mean field propagation [Eq. (||)] from to to to + r 
is calculated based on the Vlasov equation with the 
"Gaussian-Gaussian approximation" whose exact mean- 
ing is given bellow. 

We require that the Wigner function f k (x,t) of the 
nucleon k should satisfy the Vlasov equation 



St ~ »- ' «- »- ' »- ' 1 ' 



dh 
dp 



dfk 
dr 



c)h 
Or 



dfk 
dp 



at least approximately, for the time interval from to to 
to + t, with the initial condition f k (x,0) = g(x;X k ). A 
special notation of the time derivative (S/5t) is adopted 
for the mean field propagation in order to avoid a fu- 
ture possible confusion. Of course, it would be an easy 
numerical task to solve Eq. (jlT]) directly. However, in 
order to make the later decoherence process as simple as 
possible, we take another method by writing the Wigner 
function as the mean of the stochastic virtual phase space 
distributions of deformed Gaussian shape, 



f k {x,t)=g{x;X k (t),S k (t)) 



cftX 

g(x;X,S k (t))w k (X,t) — , 



(18) 



with 



g(x;X, S) 



8Vdet S 



exp 



1 

^^S^iXa-X^Xb-Xi) 

(19) 



ib=l 



Thus we have the virtually stochastic variables X ka (t) 
and S ka b(t) which represent the centroid and the shape, 
respectively, of the virtually stochastic distribution 
g{x, X k (t), S k (t)). The initial condition for them can be 
given by X ka (0) = X ka and S kab (0) = \8 ab at the ini- 
tial time t [Eq. The stochasticity of S ka b(t) is not 
shown explicitly in Eq. (|lj) because we will see later that 
its stochasticity is weak. 

We shall now determine the virtually stochastic time 
evolution of X ka {t) and S ka b(t). It is first noticed that 
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the Vlasov equation can be applied to each component 
of Eq. ( [l8|) as long as Eq. ( |l7| ) is linear in f k . The time 
evolution of gk(%, X k (t), S k (t)) by the Vlasov equation 



St 



dh dg dh dg 
dp dr dr dp 1 



(20) 



is characterized by the time evolution of the first and the 
second moments of the distribution 



(21) 



J x ka {t) = -Jx a g(x 1 t)—, 
S S f / \ / \ d^x 

J^S ka b(t) = jj; [ X a- X ka (t)J (x b - X kb (t) J g(x, t) —j- , 

(22) 

in which (S/St)g(x, t) is given by Eq. (p0[). The shape of 
the distribution at t + At is characterized by S kab (t) + 
(S/St)S kab (t)At. This symmetric matrix is diagonalized 
by an orthogonal matrix as 



S ka b(t) + -^S kab (t)At — X c O ac Obc- 



(23) 



The distribution is generally diffusing in some directions 
and shrinking in the other directions in the phase space. 
We will extract the component that is diffusing beyond 
the original width of the wave packet [Eq. (0)] by 

^Skab(s)) + = A ljni o ^ max ( ' Ac _ l) °^Obc, 

(24) 

which is now taken into account, not by changing the 
variable S ka b(t), but by giving a virtual Gaussian fluctu- 
ation AX ka (t) to the centroid X ka {t) satisfying 



AX 



AX ka {t) = 0, 
S 



(25a) 



ka{t)AX kb (t') = (-S k ab{t)) S(t-t'). (25b) 



The equation of motion for X ka (t) with virtual stochas- 
ticity may be written as 



d S 

J t Xka{t) = J t X ka {t) + AX ka (t). 

The equation of motion for S kab (t) is 



(26) 



d_ 

dt 



Skab{t) = -^S ka b(t) — (^-^S ka b(t)^ , (27) 



which does not contain the diffusing component beyond 
the original width because its effect has been counted as 
the fluctuation to X ka b(t). 

We emphasize again that the stochasticity has been, 
at this stage, introduced only as a numerical method to 
solve the mean field propagation [Eq. (pTj)]. By tak- 
ing the mean of the ensemble of stochastic distribu- 
tions g(x,X k (t), S k (t)) [Eq. (|l8|)l, our solution will re- 
produce the deterministic solution of Eq. (J17l), up to 



the Gaussian-Gaussian approximation. This approxi- 
mation means that we have introduced some restriction 
on f k (x,t) by considering only the Gaussian fluctua- 
tion to the centroid of deformed Gaussian distribution 
g(x, X k (t), S k (t)). Nevertheless, f k (x,t) is not restricted 
to a Gaussian form because the stochastic centroid X k (i) 
can move on its own way in each stochastic realization. 



3. Decoherence 

At the time to+r, we will now perform the decoherence 
into AMD wave functions [Eq. (jg)] , by decomposing the 
wave function of each nucleon fc, 



\Mt))(Mt)\ 
(Mt)\Mt)) 



\ip(w)) (ip(w)\ 
{tp(vf)\<p(w)) 



w' k (w,r)dw, (28) 



where |V'fe( r )) stands for the state after the mean field 
propagation from to to to + t with the initial condition of 
the Gaussian packet \ip(W k )) at to- In the phase space 
representation, an equivalent equation is written as 



fk(x,r) 



g{x;X)w k (X,r) 



d b X 



(29) 



with the wave packet g(x; X) defined by Eq. (|14j). We are 
going to show that the weight w' k (X, r) in this equation 
is given by the weight w k (X, r) defined by Eq. ( p"8| ) under 
a reasonable choice of the requirement. 

Before deriving it, we should notice two features of 
the variable S ka b(T~) which represents the shape of the 
stochastic Gaussian distribution g(x; X k (r), S k (r)). The 
first feature is that the stochasticity of S ka b{T) is weak 
because it is only through the indirect influence of the 
stochasticity of X ka (t) [Eq. ( p7|) and Eq. (gq)], and there- 
fore 

S kab (r) w S' kab (r) if \X ka (r) - X' ka (r)\ < 1, (30) 

for different stochastic realizations (X ka (t), S ka b{t)) and 
C^fca(*)'^ka&(*))- The second feature is that the eigen- 
values of S kab (T) (Ai < • • • < Ae) are bound by and | 
as is evident from the way of the construction [Eq. fl27Q 
with Eq. (p4j)] and normally three of them are equal to 

4" 



< Ai < A 2 < A 3 < X 4 



A s 



A 6 = 2 . 



(31) 



[These two features are valid if the coherence time r is 
not very long or the mean field Hamiltonian is not very 
different from a quadratic form (with arbitrary curva- 
tures) in the interesting phase space region, as we will 
see later in simple examples.] For convenience, we can 
change the variable from x to y, the latter being defined 
so as to diagonalize S kab (r), 



Vb = y^Oqfc^a, 



Y b — ^ O ab X ai 



(32) 
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where O a b is the matrix to diagonalize Skab(j)i 

Skab(r) = ^ ^cPacObc- (33) 

c 

We now require that, integrated out the shrinking 
directions, the distribution in the diffusing directions 
(2/4,?/5, 2/6 ) 

dy 1 dy 2 dy 3 



fk(y,r) 



r 3/2 



(34) 



should be unchanged by decoherence. From Eqs. ( |18| ) and 
(p9|), it is easily seen that the requirement is satisfied by 
choosing w' k (Y. r) so that 



w' k (Y,r) 



dY 1 dY 2 dY 3 

^3/2 



w k (Y, 



dY 1 dY 2 dY 3 



(35) 



This condition is, of course, fulfilled by taking w'(Y, r) — 
w(Y,t). 

What about for the shrinking directions (jji, y 2l 2/3)? In 
typical cases, the width of the Wigner function fk{y, t) in 
the shrinking directions is comparable to Ai,A2,A3 and 
smaller than j . [This is because the fluctuation has been 
given only to the diffusing directions and therefore the 
weight distribution w(Y, t) is narrow in the shrinking di- 
rections, which is the case if r < 2n/ui, where us is the 
oscillation frequency corresponding to the curvature of 
the mean field Hamiltonian.] Therefore, by using Gaus- 
sian packets of width | and positive weight w' k in Eq. 
(p9|), it is impossible to reproduce the shrinking compo- 
nent of /fc(2/,T). This is a physical consequence of the 
decoherence into wave packets. After decoherence has 
taken place, the shrinking disappears due to the uncer- 
tainty principle in each branch, which should be regarded 
as a physical change of the state by decoherence. There- 
fore we shall never require that the distribution in the 
shrinking directions (2/1,2/2, 2/3) would be kept unchanged 
by decoherence. Instead, we require that the width of 
fk{y,T~) in these directions is kept as close as possible to 
the minimum value j. As mentioned above, the weight 
distribution w(Y, r) in the shrinking directions is usually 
narrow, and therefore the choice w' k (Y, r) = w(Y, r) does 
not increase the width of fk(y,i~) in the shrinking direc- 
tions much beyond |. 

The derived numerical procedure for decoherence is 
quite simple. The virtually stochastic variable Xfe a (r), 
obtained by the mean field propagation procedure, is now 
given the physical meaning as the wave packet centroid. 
Accordingly, the variable of the shape is replaced as 



Skab(r) — > J^ab, 



(36) 



and the calculation is continued to the next mean field 
propagation. 



depends on the state of the other nucleons. If one follows 
the original idea of the mean field propagation, h should 
be calculated as h[f(t)] by using the Wigner function 
fi(x,t) given by equations similar to Eq. ( |l8|) for I = 
1,. .. ,A. At the time to + t when decoherence takes 
place for single particle wave functions, the mean field 
Hamiltonian will also change as 



M/»] - h\g(X(r))], 



(37) 



the latter being calculated for the wave packets 
g(x, X[(t)). We will call this change "mean field branch- 
ing" . 

We may think of the other possibility that mean field 
branching does not necessarily take place at the same 
time as the decoherence of single particle wave functions. 
Let us introduce the time scale r m f of mean field branch- 
ing which can be different from r, and generalize Eq. (|3T 
to 



hit) I WW 
ft W " i h[g(X(t))] 



for 
for 



< t < T mf 
Tmf <t <T 



(38) 



In fact, in a more elaborate theory, the decoherence of 
wave functions may take place gradually during the time 
interval from to to to + r. Then it can be reasonable 
in our model with sudden branching to take r m f shorter 
than r; for example, r m f = ^t. As mentioned in Sec. II A, 
another physical possibility of r m f < r is found if one con- 
siders the decomposition like Eq. (||) before the coherence 
time r. Even though the pre-branches have not deco- 
hered, the mean field approximation may be applicable 
not to the total state 1^) but to each of the pre-branches 
|<//) (g) |#'), \tp") <g) |#"), . . . separately. This situation 
is possible because the true time evolution is linear while 
the mean field approximation is nonlinear. If this is the 
physical case, we should take r m f < r. It should also 
be mentioned that, if the coherence time r is short com- 
pared to the time scale of the diffusion l/[(5/5t)Si] + , the 
result does not depend on the choice of T m f (< r). 

The extreme case of r m f = is convenient for the 
numerical calculation with a code based on molecular 
dynamics. In this case, we only need to calculate the 
mean field Hamiltonian h g(X(t))] without performing 
the mean averaging in Eq. (lq) which would be a hard nu- 
merical task. Another merit is that h[g(X(t))] can be re- 
placed with a precise mean field Hamiltonian /i[$(Z(i))] 
which is obtained from the fully antisymmetrized AMD 
wave function \(fr(Z(t))) without employing the physical 
coordinate approximation. 



5. Equation of motion and energy conservation 



4- Mean field branching 

Equation (^) for the mean field propagation of the 
nucleon k contains the mean field Hamiltonian h which 



The equation of motion ( J2q ) should be written in the 
original coordinates Z . Furthermore, a special consider- 
ation is necessary so as to ensure the total energy con- 
servation. No change has been made since Ref. |10| and 
these problems are not directly related to the main aim of 
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this paper. Therefore the readers who are not interested 
in them can skip this part and go to Sec. |ll C . The pre- 
scriptions which have been given in Ref. |10(| are briefly 
summarized bellow. 

Before writing down the stochastic equation of motion 
for the wave packet centroids Z, several comments are 
necessary. The deterministic part of the equation is de- 
rived by the time-dependent variational principle for the 
AMD wave function rather than using the deterministic 
part of Eq. ( |26| ) , though these two ways should be almost 
equivalent. The fluctuation AXk a is for the physical co- 
ordinate Xka and therefore it is necessary to convert it 
to the fluctuation for the original coordinate Z. This is 
done by introducing a stochastic one-body quantity Ok 
that generates the fluctuation AXk in the form of the 
Poisson brackets, as shown in Ref. |ic] . It has been im- 
plicitly assumed that there is no correlation among the 
fluctuations of different wave packets. However, some 
minimal correlations should exist so that the conserva- 
tion laws are satisfied. 

The equation of motion for the centroids is written as 



d r 

m k 

PkrnQ III. 



k=l 



(39) 



where the Poisson brackets {J- ', G} and the inner product 
of the gradients (J 7 , Q) are defined by 

dg dg x dT 
ih ^\dz^^ T dz* T dz ia ; ^dz* T 



(40) 



<&_ c -i dg , dg „_ x dT 



1(TJT 



dZ* 



+ 



dZi 



-CI 



dZ* 



(41) 



dz* a dz JT 



log($(Z)|$(Z)) a,T = x,y,z. 



(42) 



The subscripts Ck and Nfc attached to these brackets 
indicates that the consideration is limited to the sets of 
nucleons Ck and respectively, where Cfe stands for 
the cluster that includes the nucleon k and stands for 
a neighborhood of the nucleon k. The explicit definition 
is given in Ref. Q . 

The first term of Eq. ([59|) has been derived based on 
the time-dependent variational principle, with TL given 
by 

n(Z) = - + T (A - N F (Z)). 



2M 



(43) 



The quantum Hamiltonian H includes an effective two- 
body interaction such as the Gogny force [E5| which can 



be density dependent. The spurious kinetic energies of 
the zero-point oscillation of the center-of-mass of the iso- 
lated fragments and nucleons have been subtracted in Eq. 
( [f3| ) by introducing a continuous number of fragments 
Np(Z) J!, M. Without this subtraction, the Q-values 
for nucleon emissions and fragmentations would not be 
reproduced. The parameter To is 3h 2 is/2M in principle 
but treated as a free parameter for the adjustment of the 
binding energies. 

The first term in the square brackets of Eq. ( |39| ) is 
the fluctuation due to AXk generated by the stochastic 
one-body quantity Ok as mentioned above, with the cor- 
rection for the conservation of the three components of 
the center-of-mass coordinate and those of the total mo- 
mentum (denoted by {V m }). The Lagrange multipliers 
{a m } should be determined by 



{VuOk}c k +Y{V h V m }c k ak m = 0. 



(44) 



The second term in the square brackets of Eq. ( p^ ) is 
the dissipation term to achieve the energy conservation. 
Since the dissipation term should not violate the other 
conservation laws, the center-of-mass coordinate, the to- 
tal momentum and the total orbital angular momentum 
(denoted by {Q m }) are kept constant by determining the 
Lagrange multipliers (3k m by 



{Qi,H) Nk + V(Q/,Q m ) Nfc /?fc m -0. 



(45) 



The monopole and the quadrupole moments in the coor- 
dinate and momentum spaces are also included in {Q m } 
when Nfe is composed of more than 15 packets, because 
the global one-body quantities should be well described 
without the dissipation term due to the way how the fluc- 
tuation is derived. The parameter \ik is then determined 
by 



f'k 



{H, Ok + E, 



' - - ■■ - -tm a krnPm\Ck 



(46) 



in order to conserve the total energy TL. Finally, we need 
to avoid the problem that the fluctuation is finite even 
near the ground state and therefore the energy conserva- 
tion is impossible. This is done by introducing a reduc- 
tion factor 7^ in front of the square brackets of Eq. ( |39| ) 
near the ground state (see Ref. pOj). 

C. Specific models and simple examples 

Our framework includes two essential parameters r and 
r m f which represent the time scales for the decoherence 
into wave packets and the mean field branching, respec- 
tively. Each of the following models can be regarded as 
corresponding to a specific choice of (r, T m f) [with an ad- 
ditional simplification in the case of the original AMD]. 
One of the main aims here is to compare the last two 
models (AMD/D and AMD/DS) by taking simple exam- 
ples. 
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1. Mean field models 

It is needless to mention that the choice of r = r m f = 
oo corresponds to the mean field theory which solves the 
mean field equation [Eq. (|])] for an given initial Slater 
determinant toward the final state for large t. In fact, in 
this case, our model is equivalent to solving the Vlasov 
equation 



df _ dh[f] df | dhjf] df 
dt dp dr dr 9p' 



^ (47) 



with t he Ga ussian-Gaussian approximation introduced 
in Sec. II B 2 . The collision term || [| is not explicitly 
shown here for the brevity of presentation. 

It is worthwhile to mention the similarity and the dif- 
ference between our scheme and the mean field trans- 
port theory with fluctuation |[. In the latter theory, 
a stochastic term A/ is added to the mean field equa- 
tion. However, from the view point of this paper, we can 
equivalently interpret that the one-body Wie ner function 
/ follows the usual mean field equation (J47]) and, at the 
same time, the decomposition is made by 



/(r,p,i) = £>(A/,i) /(r,p,t) 

A/ 



A/(r,p) (48) 



which defines a scheme of mean field branching. Never- 
theless, we encounter a problem if we want to interpret 
each stochastic realization / + A/ as a decohered state. 
A special implementation of the fluctuation A/ is neces- 
sary, because the randomness of A/ does not generally 
ensure the idempotency of / + A/ and the eventual lo- 
calization of it in phase space. 



2. The original AMD 

In the original version of AMD H, the change of the 
wave packet shape in the mean field propagation [Eq. (|^)] 
has not been considered. This corresponds to replacing 
Eq. (p2|) with 



T-S'fcafc(^) — 0, 

dt 



(49) 



to have a constant shape Skab = \&ab- Then there is no 
branching due to decoherence, and we have a determin- 
istic equation of motion 



dt 



Z, — {Z i} H} 



(50) 



instead of Eq. (39). It should be noted, however, that 
the two-body collision effect has been incorporated as a 
stochastic process [|[ [l6| in addition to Eq. ( |50| ) already 
in this original version, as well as in all the other versions 
of AMD. 




FIG. 1: The branching of the wave packet in AMD/D is 
schematically shown for a free nucleon. The <-» symbols show 
the fluctuation to the wave packet centroids. Light gray region 
shows the exact time evolution of the Wigner function /. 



3. AMD/D 

In Refs. §, [n}, the wave packet diffusion effect by 
the mean field propagation has been incorporated into 
AMD as a source of the stochastic branching of the wave 
packets. This version of AMD corresponds to the case of 
the strongest decoherence (t = r m f — * 0) in the present 
general framework. It is straightforward to confirm that 
taking the limit r — > exactly results in the formulation 
ofRef. |TJ. 

This version of AMD, which is also called AMD-V, 
is called AMD/D in this paper, because the wave packet 
diffusion (D) effect in the mean field propagation [Eq. (^) 
or Eq. (O)] is taken into account as branching while the 
shrinking of wave packets in some phase space directions 
is discarded. It should be noted that the shrinking is 
a result of the coherent mean field propagation which 
disappears in the case of the zero coherence time r = 0. 

It is difficult to judge a priori whether the choice of 
t = is reasonable or not. The answer will depend on 
the type of reaction which we want to describe. AMD/D 
has been applied, with a reasonable success, to various 
reaction systems in the medium energy region not only 
for nuclear collisions ||, [TTJ, |2tj, pljbut also for nucleon in- 
duced fragmentation reactions |[22|. However, these suc- 
cesses do not mean that AMD/D is valid for all kinds of 
reactions. 

We can a priori expect that AMD/D will fail if deco- 
herence is not a physical case. The most simple and clear 
example is the dynamics of a single nucleon moving in a 
one-body external potential U (r) , in which decoherence 
does not exist because the nucleon is not interacting with 
anything. In this case, the exact time-evolution ^(r, t) is 
given by the time-dependent Schrodinger equation with 
an initial condition ip(r, to) cx exp[— v(r — Z/y^) 2 ]. It is 
known that Eq. (|TJ) with h{r,p) = p 2 /2M + U(r) yields 
the exact quantum-mechanical time evolution of the one- 
body Wigner function /(r, p, t) when the potential has a 
quadratic form, including the case of a free nucleon or a 
nucleon in a harmonic oscillator potential with arbitrary 
curvatures. 

The light gray region in Fig. [l] shows the exact time 
evolution of the Wigner function /(r, p, t) for a free nu- 
cleon with the initial condition of a Gaussian wave packet 
at to. The spatial distribution increases as the time 
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progresses, while the momentum distribution does not 
change at all. Due to the Liouville theorem, the phase 
space volume is conserved, which is the semiclassical ana- 
logue of the fact that /(r, p, t) corresponds to a pure state 
ip(r, t). At the initial time to, the Wigner function is dif- 
fusing in the direction of the <-> symbol in Fig. [I] and 
shrinking in the other direction. If AMD/D is applied 
to this situation, the diffusion is taken into account as 
branching by giving fluctuation to the wave packet cen- 
troid in the «-> direction, while the shrinking is ignored. 
At another time, t\ or i 2 > each of the branched wave pack- 
ets is treated in the same way as the initial wave packet 
at to (except for the different centroid value), without 
any influence of the history of the wave packet and the 
existence of the other branches. Namely, the fluctuation 
always has the same property for an free nucleon, and 
therefore both the spatial and the momentum distribu- 
tions increase as the time progresses. Consequently the 
coherent Wigner function /(r, p, t) cannot be reproduced 
if AMD/D is applied. Generally speaking, decoherence 
increases the width of the phase space distribution. 

In the practical calculations of nuclear reactions, the 
above problem for a free nucleon is not so serious because 
the dynamics of free nucleons is not of our interest, and 
it is usually sufficient to describe an emitted nucleon as 
moving on a straight line (or on a Coulomb trajectory) 
as a classical particle. Therefore, in Refs. [[)], [l0| the 
branching was switched off for isolated nucleons. For 
example, the switching-off condition adopted in Ref. |l(| 
is 



to 



1.75 



(Zi-Zfc)l < 10 



(51a) 



and 



1.75-|R(Zi-Z fc )| R(Zi-Zfc) 



< 5. (51b) 



This condition switches off the branching for the nucle- 
ons in small clusters with A c < 10 as well as emitted 
nucleons. 

Not only the branching is switched off for isolated nu- 
cleons, but also each isolated nucleon is regarded as hav- 
ing a definite momentum value without the internal dis- 
tribution of the Gaussian wave packet. This interpre- 
tation is consistent to the definition of the Hamiltonian 
[Eq. (ff3j)] where the zero-point oscillation kinetic ener- 
gies of isolated nucleons have been subtracted. Therefore, 
we can get reasonable results even with AMD/D if the 
switching-off condition is appropriately chosen so that 
the branching is switched off (at t 2 in Fig. [j] for example) 
when the momentum centroid has got the appropriate 
amount of the fluctuation corresponding to the internal 
momentum distribution of the initial wave packet. 

However, it will not be easy to find the switching-off 
condition that works for all situations. The condition 
of Eq. (^TJ) may not work well for the reaction systems 
which have not been studied. When a big system is ex- 
panding slowly, for example, the switching-off will not 




FIG. 2: The branching of the wave packet in AMD/DS 
is schematically shown for a free nucleon. Light gray re- 
gion shows the exact time evolution of the Wigner func- 
tion /. Each of the solid ellipses shows the shrunken shape 
g(r, p; X(t), S(t)) in each branch, while the dashed circles 
show the wave packets with the original width. The <-» sym- 
bols schematically show the magnitude and the direction of 
the fluctuation to the wave packet centroids. 



take place for a long time, and then the branching will 
continue too long for the AMD/D method to reproduce 
the mean field prediction. In fact, we will show in Sec. 
Ill that the AMD/D method with the switching off con- 
dition (|5l|) seems to overestimate the diffusion effect in a 
rather slowly expanding big system. Therefore we want 
such a new scheme of branching that we have no ambi- 
guity of introducing any switching-off condition. 



4. AMD/DS 

Now we introduce a new scheme of decoherence by tak- 
ing the choice of a large coherence time r and a short time 
scale for mean field branching r m f = 0. We will call this 
model AMD/DS because the wave packet shrinking (S) 
effect in the mean field propagation [Eq. ([|) or Eq. P?])] 
is reflected in the dynamics as well as the diffusion (D) 
effect with the choice of a finite t. The explicit defini- 
tion of r, which depends on two- nucleon collisions, will 
be given later. 

AMD/DS provides us a one-body dynamics similar to 
that in mean field models, as expected from the fact that 
the coherent mean field propagation is respected by the 
choice of a large coherence tim e r. Figure || illustrates 
how our formulation in Sec. II B works in the simplest ex- 



ample of a free nucleon. The light gray region is identical 
to that in Fig. [|, showing the exact time evolution of the 
Wigner function /(r, p, t) with the initial condition of a 
Gaussian packet at the time to. The wave packet diffu- 
sion effect, [(S/St)S(t)] + defined by Eq. (gj), is taken into 
account as before by giving the fluctuation to the wave 
packet centroid X(t) = (R(t),P(t)) in Eq. @ with the 
property of the fluctuation given by Eq. (|25f). On the 
other hand, the shrinking effect is reflected to the shape 
matrix S(t) by the equation of motion (|27|). The ellipse 
with a solid line in each branch in Fig. ^represents the 
deformed and shrunken shape of g(r,p; X(t), S(t)). It 
should be noted that the eigenvalues of S(t) do not ex- 
ceed the original wave packet width \ (shown by dotted 
circles in Fig. @) since the diffusion beyond that width is 
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considered by the fluctuation to the centroid. The exact 
Wigner function /(r, p,t) is reproduced by the mean av- 
erage of the elliptic shapes g(r,p;X(t),S(t)) as defined 
by Eq. £§. 

A general difference between AMD/DS and AMD/D 
is that the phase space diffusion is weaker in AMD/DS 
than in AMD/D, reflecting the different strength of de- 
coherence. Mathematically, this difference arises due to 
the property of the fluctuation to the wave packet cen- 
troids. Namely, the Vlasov equation (|2(]) is applied to 
different g(r, p; X(t), S(t)) which has always a full width 
in AMD/D and has a shrunken shape in AMD/DS, and 
therefore the property of the fluctuation [Eq. (pq)] is dif- 
ferent. In general, AMD/DS has smaller strength of the 
fluctuation than AMD/D (except for the switching-off in 
AMD/D). In the case of the AMD/DS description of a 
free nucleon in Fig. g, as the time progresses, the strength 
of the fluctuation gets smaller and smaller with the di- 
rection of the fluctuation also changing, and at t = oo 
the fluctuation ceases and the elliptic distribution iS(oo) 
is completely shrunken to have a definite value of the 
momentum. 

It is easily proved that, in some cases, the coher- 
ent time evolution of the Wigner function is repro- 
duced exactly by AMD/DS as the mean average of the 
shrunken elliptic shapes g(r, p, X(t), S(t)), in spite of the 
Gaussian-Gaussian approximation adopted in Sec. II B 1 . 



This is the case not only for a free nucleon but also for a 
nucleon in a harmonic oscillator potential with arbitrary 
curvatures. In general cases, however, the reproduction is 
not exact because of the Gaussian-Gaussian approxima- 
tion. In future, such an approximation may be removed 
if it turns out to be necessary, but in the present work 
we are satisfied with this agreement of AMD/DS, which 
ensures the approximate validity of AMD/DS in the ex- 
treme case when decoherence does not physically take 
place. 

When the coherence time r is long and the mean field 
branching time is short (r m f = 0), it is interesting to 
imagine the behavior of this model in nuclear collisions. 
We can expect that the mean field branching (namely, 
the stochastic fluctuation of the mean field) does not in- 
fluence so much the global dynamics and the early dy- 
namics in high density stage. For such aspects, AMD/DS 
will behave similarly to the mean field model because the 
coherence of the single particle wave functions is kept for 
a large time scale r. On the other hand, for the aspect 
of clusterization, AMD/DS with r m f = will behave like 
usual molecular dynamics models, because the branched 
mean field h[$>(Z(t))] is equivalent to the mean field in 
such models, and therefore it helps the formation of clus- 
ters each of which is bound by the branched mean field. 

The policy of the scheme of AMD/DS is to take r that 
is as large as possible. However, if a two-nucleon collision 
takes place, it will not make sense to keep the coherence 
of the single particle wave functions of collided nucleons. 
Therefore, we assume that decoherence takes place for 
a nucleon with some probability when it experiences a 



two-nucleon collision with another nucleon. [Therefore r, 
defined for each nucleon, is the time interval between two 
successive collisions related to it.] The probability Pdec 
of decoherence at each two-nucleon collision is chosen to 
be 



Pdec(E,t 



,-E(l-coa6)/Eo 



(52) 



where E is the two-nucleon collision energy in the labo- 
ratory system for the two nucleons and 9 is the scattering 
angle in the center-of-mass system for the two nucleons. 
The purpose of this probability is to reject the low mo- 
mentum transfer cases where the scattered state has a 
significant overlap probability with the case of no col- 
lision. Note that the probability is related to the mo- 
mentum transfer q 2 = ME(1 — cos 8). The parameter 
Eq = 15 MeV is chosen in the present work. With this 
choice, decoherence takes place in most of the collisions 
between the nucleons from the different initial nuclei in 
the early stage of the reaction with the incident energy 
more than several ten MeV/nucleon, while decoherence 
seldom takes place within the initial nuclei and the pro- 
duced clusters. 

At the end, let us discuss on the method of the sub- 
traction of the zero-point kinetic energies of the isolated 
nucleons in Eq. ([l3|). This subtraction is consistent to 
the coherent one-body dynamics of AMD/DS, in that the 
shrunken shape g(r, p; X(t), S(t)) for an emitted nucleon 
eventually has a definite momentum and therefore the 
zero-point kinetic energy should not be counted in the 
conserved energy Ti. However, when a nucleon is coming 
out of a nucleus, the zero-point energy subtraction plays 
as a repulsive force to the nucleon from the nucleus in 
spite of the fact that such a repulsion does not exist in 
the one-body dynamics in the mean field. It is possible to 
remove the repulsive effect by keeping the conserved en- 
ergy Ti. given by Eq. ( |43| ) . This is easily done by formally 
adding a term to the fluctuation part in the equation of 
motion (B9h by the replacement 



dZt 



dO k dN F 
+ Oj k T - 



dZ. 



dZ* 



(53) 



Then the added term cancels the zero-point subtraction 
term in the deterministic part {Zi,7i} and thus the re- 
pulsive force mentioned above does not exist any longer. 
The replacement of Eq. ( p3|) is done at any place where 
Ok appears, like in Eqs. (|44|) and (46). Therefore, the 
conserved energy is still Ti with the zero-point energy 
subtracted. The key of this trick is that the zero-point 
energy is converted to the translational motion of the nu- 
cleon in the old treatment while that energy is shared by 
all the nucleons in the new treatment with the replace- 
ment of Eq. (53). 
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III. EFFECTS IN A MULTIFRAGMENTATION 
REACTION 

In this section, we discuss the results for 129 Xe + Sn 
collisions at the incident energy E/A = 50 MeV and 
the impact parameter range < b < 4 fm. A detailed 
and systematic analysis of this reaction system will be 
given in separate papers p3| , p4| . The main purpose here 
is to make comparison of the two models of the quan- 
tum branching in AMD, one (AMD/D) with only the 
wave packet diffusion effect and the other (AMD/DS) 
with the wave packet diffusion and shrinking effects, the 
latter effect being a consequence of the coherent mean 
field propagation for a finite coherence time r. From the 
character of these models, as discussed in the previous 
section, we expect that differences should be found in 
the diffusion property of nuclcons in nuclear matter and 
the global one-body dynamics. 

Many events with various impact parameters in the 
range of < b < 4 fm were produced by solvi ng the 
stochastic equation of motion given in Sec. II B. The 
triple loop approximation flfjfl was used in order to save 
the computation time. The Gogny force |2^] was used as 
the effective interaction. In the calculation of AMD/DS, 
we use the new treatment of the subtraction of the zero- 
point kinetic energies of nucleons and clusters given at 
the end of Sec. II C 4, while AMD/D calculation is done 



with the old treatment of the zero-point subtraction, un- 
less otherwise stated. In addition to the equation of mo- 
tion, the two-nuclcon collision effect was introduced in 
the usual stochastic way j^, The two-nucleon colli- 
sion cross section adopted here is given by 



cr{E,p) = min a LM (E,p) 



100 mb 



1 + £/(200 MeV) 



(54) 



where ctlmCE;, p) is the cross section given by Li and 
Machleidt [M from Dirac-Brueckner calculations using 
the Bonn nucleon-nucleon potential. This cross section 
depends on the two-nucleon collision energy E and the 
density around the collision point p. It also depends on 
the isospins of the colliding nucleons. The temperature 
in the parameterization by Li and Macheit was, however, 
replaced by zero. A low energy cut has been introduced 
in the adopted cross section, though its effect on the final 
results has turned out to be unimportant. For the angu- 
lar distribution, we use the same parameterization as in 
Ref. jl6). 

The calculation of each event was started by putting 
two nuclei with a distance 12 fm and boosting them at 
the time t — 0. The AMD calculation was continued 
until t — 300 fm/c, at which we assume the thermal 
equilibrium of each produced fragment and calculate its 
decay by using a statistical decay code Q which is based 
on the sequential binary decay model by Piihlhofer f p7| . 

Figure |3J shows the time evolution of a typical event 
with the impact parameter b = 3.4 fm. It appears that a 
system is formed, which, after a maximum compression 
around t ~ 45 fm/c, expands and many clusters appear 



around t ~ 100 fm/c. The expansion is stronger in the 
beam direction than in the transverse directions, which 
means that the initial nuclei do not stop completely even 
in such central collisions. Therefore, another possible in- 
terpretation may be that the initial two nuclei are pass- 
ing each other with large dissipation and breaking up 
into clusters. However, the mixing of the wave packets 
from the two nuclei is considerable. On the average, 87 
nucleons from the projectile 129 Xe nucleus come out to 
the forward direction (p z > in the center-of-mass sys- 
tem), while the other 42 nucleons from the projectile nu- 
cleus appear in the backward direction. This corresponds 
to around 67%-33% sharing of the projectile nucleons in 
forward-backward directions, to be compared with 50%- 
50% for full mixing or 100%-0% for no mixing at all. 
These qualitative features do not depend so much on the 
choice of the models of the quantum branching, though 
the event of Fig. |^ was obtained with AMD/DS. 

The same reaction system has been studied by Nebauer 
et al. with the quantum molecular dynamics (QMD) p8[ . 
A serious problem of their QMD result is that too large 
projectilelike and targetlike fragments are produced even 
in the central collisions (i?trans > 450 MeV). Conse- 
quently the QMD calculation largely overestimates the 
yield of the big clusters with Z > 20, as shown in Fig. 7 
of Ref. [^8| . This problem of the spurious binary feature 
is qualitatively similar to the problem which has ever en- 
countered in the AMD calculation by Ono and Horiuchi 
for the 40 Ca + 40 Ca collisions at 35 MeV/nuclcon [§. 
This problem in AMD has been solved in Ref. (|] by the 
stochastic incorporation the wave packet diffusion effect 
which allows the mixing and/or breakup of the initial nu- 
clei. In fact, the present AMD calculation do not show 
the binary feature as strong as in the QMD calculation, 
which can be seen in Fig. ^ and in the cluster charge dis- 
tribution to be shown later. The fermionic nature may 
also be important to solve the problem of QMD. Papa et 
al. have shown in Ref. p9| that the spurious binary fea- 
ture disappears when they introduce the fermionic nature 
into QMD in a stochastic way. 

We can expect that the early stage dynamics is not 
so sensitive to the model of quantum branching, be- 
cause AMD/D and AMD/DS are equivalent for a short 
time scale and because the effective coherence time in 
AMD/DS is short due to many two-nucleon collisions. 
In fact, it is seen in the part of t < 60 fm/c in Fig. |J, 
which shows the time dependence of the two quantities 



? 2 +R 2 

Hx ' iy 



.1 



(55) 



(56) 



characterizing the transverse expansion. We use the 
transverse components of the physical positions R; and 
the physical momenta defined by Eq. (|l2|) . The brack- 
ets stand for the averaging over the events with the im- 
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FIG. 3: The time evolution of the density in the center of mass system projected on to the reaction plane, in a typical event 
of the 129 Xe + Sn collision at 50 MeV/nucleon, from t = to t = 225 fm/c. The beam direction is parallel to the horizontal 
axis, and the impact parameter of this event is 3.4 fm. The size of the shown area is 60 fm x 60 fm. This is a result of the 
calculation with AMD/DS. 
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FIG. 4: Time evolution of the transverse radius {R±) (solid 
lines) and the transverse momentum (P±) (dashed lines) in 
129 Xe+ Sn collisions at 50 MeV/nucleon averaged for the im- 
pact parameter region < b < 4 fm. The results of AMD/D 
are shown by thin lines, while those of AMD/DS are shown 
by thick lines. Filled and open circles show (R±) and (Px), 
respectively, obtained by the origial AMD calculation applied 
to the intermediate states at t = 60 fm/c of the AMD/DS 
calculation. 



pact parameter range < b < 4 fm. [We consider these 
quantities rather than the root mean square quantities 
in order to focus on the central part of the system where 
clusters are mainly produced.] The solid lines in Fig. 
|] show the transverse radius (R±) and the dashed lines 
show the average transverse momentum (-Pi). The thin 
lines show the result of AMD/D, while the thick lines 
show the result of AMD/DS. We can see that the trans- 
verse momentum is produced in the early stage of the 
collision (before 40 or 50 fm/c). As expected, there is no 
significant difference between AMD/D and AMD/DS for 
the early stage dynamics t < 60 fm/c. 

Now our interest is in the evolution of the expanding 



system which has been created by the early stage dy- 
namics before t ~ 60 fm/c. In Fig. ||, we can see that 
an important deviation between the two models appears 
in the spatial radius (R±) in later reaction stage. The 
expansion velocity is slower in AMD/DS with the wave 
packet shrinking effect than in AMD/D without it. On 
the other hand, the transverse momentum, which is al- 
most constant for t > 60 fm/c, is almost independent of 
the model of the quantum branching. It should be noted 
that the expansion is governed not only by the momen- 
tum centroids of wave packets but also by the property of 
the fluctuation to the wave packet centroids. The latter 
is the difference between the models, which can be natu- 
rally understood because the wave packet shrinking effect 
reduces the strength of the fluct uation to the wave packet 
centroids, as discussed in Sec. II C 4. To see this effect 



more clearly, we also show, by filled and open circles in 
Fig. ^, the results of the original AMD without quantum 
branching. In order to avoid the influence of the different 
early stage dynamics, the original AMD was applied to 
the intermediate states at t = 60 fm/c of the AMD/DS 
calculation. We can see that the transverse expansion is 
very weak without quantum branching. 

In order to get a deeper understanding, let us first con- 
sider how the expansion dynamics is described if a mean 
field model (such as TDHF) is applied. If the two-nucleon 
collision effect is negligible in the expanding system, most 
of the single particle wave functions will widely spread 
over the space. Clusters will not be produced in a mean 
field model but the global one-body distribution may be 
reliable. Due to the coherence of the mean field propa- 
gation, the nucleon position and the nucleon momentum 
are strongly correlated in the expanded system, in a sim- 
ilar way to the case considered in Fig. || for example. If 
we focus on a local part of the expanded system, each 
nucleon has a rather sharp momentum distribution like a 
classical particle. The main aim of AMD/DS is to have 
the same global one-body distribution as in mean field 
models, when averaged over the branches. The essen- 
tial difference is that the mean field varies from branch 
to branch in the case of AMD, which is the reason why 
clusters are produced in AMD, though this difference will 
not affect so much on the global one-body distribution. 
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When clusters are formed, however, each nucleoli is local- 
ized in one of the clusters, and then it should have some 
momentum distribution that satisfies the uncertainty re- 
lation and the Pauli principle. If this momentum dis- 
tribution is considered, the global momentum distribu- 
tion will become wider than the mean field prediction. 
Now the question is whether this widening of the mo- 
mentum distribution should be respected or the coherent 
mean field propagation should be respected. When we 
apply AMD/DS, the coherent mean field propagation is 
respected and the widening of the momentum distribu- 
tion is not considered as far as the one-body dynamics is 
concerned. It may be possible that the wave packet lo- 
calization does not change the global one-body dynamics 
through complicated many-body correlations which are 
out of the scope of the present models. On the other 
hand, when we take AMD/D, we respect the widening 
of the momentum distribution due to the localization of 
the wave packet as physical decoherence, which will in- 
crease the future expansion velocity. It is not possible 
to say a priori that one model is superior to the other. 
What we can say is that AMD/DS reproduces the mean 
field prediction more precisely then AMD/D. It is an- 
other problem not discussed here whether the mean field 
prediction is always reliable or not. 

How does the different expansion velocity in the two 
models appear in the observables? As we have seen in 
Fig. |], the difference does not appear in the global trans- 
verse momentum, and therefore the energy spectra are 
not good quantities to sec the effect directly The differ- 
ent expansion velocity is not due to the different momen- 
tum but due to the different strength of the spatial com- 
ponent of the fluctuation to the wave packet centroids. 
Therefore, we should look at the quantities which carry 
the information of the increase rate of the spatial radius. 
The cluster size distribution is one of such quantities be- 
cause each cluster is formed by the nucleons with similar 
spatial positions and velocities. It has been shown that 
the cluster size decreases as the expansion velocity in- 
creases H^, [3l] . Figure || shows how the total 104 protons 
in the system are divided into clusters at t — 300 fm/c 
before calculating the statistical decay of excited clusters. 
The results of AMD/DS and AMD/D are shown in (a) 
and (b), respectively. It is clearly seen that heavy clus- 
ters are produced more abundantly in AMD/DS than in 
AMD/D, reflecting the different expansion velocity. In 
the case of the original AMD calculation circles) linked 
to the early stage dynamics of AMD/DS [shown in (c)], 
the produced clusters are much bigger than in AMD/DS, 
reflecting the very slow expansion. In Fig. ||, the final 
charge distribution after statistical decay is shown to- 
gether with the INDRA data (bullets). The data and 
the calculated results can be directly compared, since 
the filter has been applied to the calculated events in 
order to take account of the properties of the detector 
system. The result of AMD/D has a serious problem 
that the multiplicity of the heavy clusters with Z > 15 
are underestimated. Instead of the heavy clusters, the 




(a) AMD/DS (b) AMD/D 




(c) AMD/DS+AMD (d) AMD/D (4<b<7 fm) 



FIG. 5: The partitioning of the total charge into the clusters 
at t = 300 fm/c (before calculating the statistical decay of 
excited clusters) in 129 Xe + Sn collisions at 50 MeV/nucleon 
with the impact parameter < b < 4 fm. The area of each 
sector represents ZM(Z) summed over the specified region of 
Z, where M(Z) is the multiplicity of clusters with charge Z. 
(a) The result of AMD/DS. (b) The result of AMD/D. (c) 
The result of the original AMD applied to the intermediate 
states at 60 fm/c of the AMD/DS calculation, (d) The result 
of AMD/D for the impact parameter region 4 < b < 7 fm. 

relatively light clusters with Z ~ 5 ; 6 are produced too 
abundantly. Therefore, it seems that the expansion is too 
fast in AMD/D. On the other hand, the reproduction by 
AMD /DS is quite satisfactory for the charge distribution 
of the clusters with Z > 3. These results suggest that it 
is reasonable to respect the mean field prediction of the 
expansion dynamics in this reaction system which con- 
sists of rather many nucleons and is expanding with a 
moderate velocity. However, AMD/DS has a problem of 
the overestimation of the proton multiplicity and the un- 
derestimation the a particle multiplicity. Probably this 
is a side effect of the fact that AMD/DS uses the mean 
field equation so faithfully that the light cluster emission 
is not respected compared to the nucleon emission. A 
special care |}2| will be necessary to explain the direct 
production of light clusters which have only one quantum 
bound state. The use of the semiclassical version of the 
mean field equation (the Vlasov equation) may also be 
one of the reasons why the nucleon emission is overesti- 
mated. 

The above quantitative results can be affected, in prin- 
ciple, by the centrality selection method. In our calcu- 
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FIG. 6: The charge distribution of the produced clusters 
in 129 Xe + Sn collisions at 50 MeV/nucleon with the impact 
parameter < b < 4 fm, after calculating the statistical decay 
of excited clusters and applying the experimental filter for the 
detector setup. Solid histogram shows the result of AMD/DS, 
while the dotted histogram shows the result of AMD/D. The 
INDRA experimental data [G8I are shown by bullets. 



FIG. 7: Similar to Fig. ^[ Solid and dotted histograms 
are the same as in Fig. H, showing the results of AMD/DS 
and AMD/D, respectively. The stars show the result of the 
calculation with AMD/DS combined with the old treatment 
of the zero-point kinetic energy subtraction. The bullets show 
the result of the calculation with AMD /D combined with the 
new treatment of the zero-point kinetic energy subtraction. 



lation all the events with < b < 4 fm are considered 
as central events, while in experiment the central events 
are selected by using the sum of the transverse energies 
(-^trails) of observed light charged particles (Z — 1,2). 
The experimental data in Fig. || were obtained by se- 
lecting the events with the condition Etrans > 450 MeV 
f28f . [These data are identical to those which have been 
shown by the histogram in Fig. 7 of Ref. |2§| in an ar- 
bitrary scale, while in our figure they are shown in the 
absolute scale.] Nebauer et al. have shown in QMD sim- 
ulation that events with 4 < b < 6 or 7 fm are also 
mixed in the selected events with -Etrans > 450 MeV 
(Fig. 3 of Ref. Q). In order to estimate the effect of 
these semi-peripheral events, we show in Fig. ||(d) the 
charge partitioning obtained by AMD/D for 4 < b < 7 
fm. We can see that the character of clusterization in 
semi-peripheral events is not very different from the cen- 
tral events [Fig. ||(b)]. The difference between central 
and semi-peripheral events in AMD/D calculation [(b) 
and (d)] is not as big as the difference between AMD/D 
and AMD/DS in central events [(b) and (a)]. Therefore, 
a possible considerable mixture of semi-peripheral events 
in the data does not change our conclusion that too many 
small clusters are produced in AMD/D and that the re- 
production is improved by respecting the coherence mean 
field propagation as in AMD/DS. 

Finally, the dependence on the treatment of the zero- 
point kinetic energy subtraction is shown in Fig. ^ for 
the cluster charge distribution. As far as the AMD/D 
model of the quantum branching is used, the result is far 



from the experimental data irrelevantly to the treatment 
of the zero-point subtraction. In the result of AMD/DS, 
we can get a slightly better reproduction of data by using 
the new treatm ent of the zero-point subtraction given at 
the end of Sec. II C 4. It should be noted that the new 



treatment of the zero-point subtraction is more consistent 
to the philosophy of the AMD /DS model of the quantum 
branching in that the mean field prediction of the one- 
body dynamics is respected. 



IV. 



SUMMARY 



In this paper, we have given a general framework which 
determines the many-body quantum dynamics by the 
combination of the coherent mean field propagation and 
the decoherence into branched wave packets. This frame- 
work contains the mean field description and the molec- 
ular dynamics description as specific cases. The model 
given by Refs. |9[ |l(| (AMD/D) corresponds to taking 
the zero coherence time (t = 0) for the mean field prop- 
agation. In this scheme, the wave packet diffusion by the 
mean field propagation is respected by giving appropriate 
fluctuation to the wave packet centroids. However, the 
usual fluctuation was not able to describe the shrinking 
of the phase space distribution which could be respected 
only by keeping the coherence of the mean field propa- 
gation. On the other hand, we have shown in this paper 
that it is possible to implement a finite time duration 
t of the coherent mean field propagation before deco- 
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herence, even though we still adopt a branching treat- 
ment. As a consequence, in the new model (AMD/DS), 
the shrinking of the phase space distribution is respected 
as well as the diffusion. AMD/DS reproduces the ex- 
act dynamics for a free nucleon and for a nuclcon in a 
harmonic oscillator potential with arbitrary curvatures. 
In general cases, the branch-averaged one-body dynam- 
ics in AMD/DS should be much closer to the prediction 
by mean field models than in AMD/D. Nevertheless, by 
the choice of T mf = for mean field branching, clusters 
can be formed in AMD/DS as well as in usual molecular 
dynamics models, because the mean field is calculated 
with localized wave packets in each branch. The two- 
nucleon collision effect is introduced as usual, and the 
decoherence into wave packets is assumed to take place 
when a nucleon experiences a two-nucleon collision with 
a substantial momentum transfer. 

The difference of the decoherence scheme between the 
two models results in the different diffusion properties 
of nuclcons in nuclear matter and the different global 
one-body dynamics. We have applied both models of 
AMD/D and AMD/DS to the 129 Xe+ Sn collisions at 50 
MeV/nucleon in the impact parameter range < b < 4 
fm, where many clusters are produced from the expand- 
ing system with a moderate expansion velocity. The ef- 
fect of the wave packet shrinking in AMD/DS certainly 
reduces the expansion velocity compared to AMD/D. 
Reflecting this difference in the expansion velocity, the 
charge distribution of the produced clusters strongly de- 
pends on the model of decoherence into branches. With 
AMD/DS, we have larger number of heavy clusters with 
Z > 15 and smaller number of relatively small clus- 
ters with Z ~ 5,6 than with AMD/D. AMD/DS repro- 



duces the INDRA experimental data much better than 
AMD/D, which suggests that the coherent mean field 
propagation for the one-body dynamics should be re- 
spected in this reaction system where a big system is 
expanding with a moderate expansion velocity. The de- 
tailed analysis of this reaction system based on the AMD 
calculations will be given in a separate paper. 

However, we do not claim that AMD/DS is always 
superior to AMD/D or vice versa. These two models 
should be regarded as different schemes of approxima- 
tion. AMD/DS respects the coherent mean field propa- 
gation, while AMD/D respects the existence of strong 
many-body correlations which causes the decoherence 
into branched wave packets. Although the decoherence 
has been considered in AMD/DS based on the two- 
nucleon collisions in this paper, it is also possible to 
have other many-body effects to cause decoherence, with 
which AMD/DS can be closer to AMD/D depending on 
the considered reaction systems. In future works, it will 
be important to investigate such possibilities. 
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